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Abstract 



(N 
> 

^ The rare K — > itvv decays play a central role in testing the Standard Model 

Q\ and its extensions. Upcoming experiments plan to measure the decay rates with 

high accuracy. Yet, unknown higher-order electroweak corrections result in a size- 
able theory error. We remove this uncertainty by computing the full two-loop elec- 
troweak corrections to the top-quark contribution X t to the rare decays Kl — > 7t°vv, 
K + —7- it + vD, and B — > Xd s vv in the Standard Model. The remaining theoret- 
ical uncertainty related to electroweak effects is now far below 1%. Finally we 
update the branching ratios to find Bt(Kl — > tt^vu) = 2.43(39) (6) x 10 -11 and 
Br(K + — > ir + vv) = 7.81 (75)(29) x 10 -11 . The first error summarises the parametric, 
03 the second the remaining theoretical uncertainties. 



1 Introduction 

The branching ratios of the rare K + — > 7i + uu and Kl — > 7i°uu decays are dominated by 
contributions of internal top-quarks in the Standard Model. This short distance sensitivity 
results in a precise theory prediction, but also in a proportionality to powers of V t *V t d- 
Accordingly, the branching ratios are suppressed with respect to generic new physics sce- 
narios by the near diagonality of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. This 
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leads to a high sensitivity to new physics, and a precision measurement of these modes 
could provide a decisive test of the Standard Model and its extensions. 

This potential will be exploited by a new generation of experiments (NA62 at CERN, 
KOTO at JPARC, and the proposed future experiment P996 at Fermilab), which aim at 
measuring the branching ratios with unprecedented precision. 

In the Standard Model the K — > ixvv decays proceed through Z-penguin and elec- 
troweak box diagrams which exhibit a power-like GIM mechanism. This implies a suppres- 
sion of non-perturbative effects and, related to this, that the low-energy effective Hamilto- 
nian [1, 2] 



involves to an excellent approximation only a single effective operator. Here Gp is the Fermi 
constant, a the electromagnetic coupling and #w the weak mixing angle. The sum is over 
all lepton flavours, A; = V* s Vid comprise the CKM factors, and fi represents left-handed 
fermion fields. 

The functions X 1 constitute the charm-quark contribution to T-L c s and add 30% to 
the total branching ratio of the K + — > -n + vv decay, while they leave the CP violating 
Kl — > 7r°z/z/ decay unaffected. The theoretical uncertainty in X 1 is 2.5% after next-to- 
next-to-leading order (NNLO) QCD [3, 4, 5] and next-to-leading order (NLO) electroweak 
corrections [6] are taken into account, and the resulting error in the branching ratio is 
small. 

The situation is different for the function X t which includes internal top-quark loops: 
it gives either the sole or the dominant contribution to the neutral or the charged de- 
cay modes, respectively. A two-loop electroweak calculation should cancel the sizeable 
scheme dependence of the input parameters. Yet, only NLO QCD corrections [2, 7, 8] 
and the leading term of the large-mt expansion of the two-loop electroweak corrections are 
known. While unknown higher-order QCD corrections result in a 1% uncertainty in X t , 
the uncertainty related to unknown sub-leading electroweak contributions is estimated to 
be ±2% [9]. This can be understood in the following way: the matching calculation with 
internal top-quark loops is purely short distance, the resulting operator renormalises like 
a current, such that the QCD perturbation theory converges well. Yet the on-shell scheme 
counterterm of sin^w includes large higher terms in the large-m t expansion. Hence the 
renormalisation scheme dependence of a/ sin 2 9\y in (1.1) cannot cancel if only the leading 
term in the large-mt expansion is taken into account. This was found in Reference [9] where 
the scheme difference between the on-shell scheme and the MS scheme was only decreased 
from 5.6% to 3.4% through the inclusion of the first order in the large-mt expansion. 

In this paper we will improve on the analysis of Reference [9] and compute the full 
electroweak two- loop corrections to the top-quark contribution X t . Only in such a way 
is it possible to fix the definition of the electroweak input parameters and reduce the 
uncertainty due to unknown higher order electroweak corrections from 2% to the per mil 
level. Since a 2% uncertainty in X t scales up to a 3% to 4% uncertainty in the branching 




(1.1) 



2 



ratios such a reduction of the theoretical error is important in particular in light of the 
coming experiments. In addition, our results are equally applicable for the B — > X^^vP 
decays. 

Our paper is organised as follows. In Section 2 we discuss the dependence of our result 
on different renormalisation schemes. In Section 3 we present some technical details of 
our calculation. Our numerical results are contained in Section 4. In the Appendices we 
provide the analytic form of the electroweak correction to X t in the limit of small sin#v^ 
and compare our expansion for a large top-quark mass with the literature. 



2 X t beyond leading order 

The truncation of the perturbation theory results in a residual scale and scheme dependence 
of the matrix elements of the effective Hamiltonian in Equation (1.1). For the top-quark 
contribution, the matrix element of the operator 

Qu= ^LlM{y lLl ^v lL ) (2.1) 

l=e,/j,,r 

factorises and 4aGp/(2\/27rsm 2 9\y)\X t will be independent of the renormalisation pro- 
cedure after higher-order corrections are included. Let us now discuss the dependence on 
the electroweak renormalisation scheme and how to combine these schemes with the NLO 
QCD results, which are known in the MS scheme. 

Pure QCD corrections leave Gf, a, and sin 2 9w unaffected, such that X t is a renormal- 
isation scheme invariant quantity if electroweak effects are ignored. It is then customary 
to expand 

X t = X t (0) + ^X 4 (1) + ^Xf w) . (2.2) 

4ll 4:71 



in terms of the leading-order (LO) contribution [10] 

(2.3) 



X t (0) = - 



x t + 2 3x t - 6 

+ 7 7^7 ln x t 



x t -l (x t - l) 2 



where Xt = mf/Myy- The schemes for m t and Mw are defined below. The NLO QCD 
correction [2, 7, 8] 



(i) 29x t — x\ — Ax\ x t + 9xf — xl — xf 

X-i = r— — ; r-r In Xf 

1 3(1 -x t ) 2 (l-x t ) 3 

8x t + 4.x 1 . + x, 3 - xf , , 4x t - x, 3 T . , <9X, (0) , uf 



(2.4) 



2(1 -x,) 3 " l (l-xtf z y t! 1 dx t M 2 



fixes the renormalisation scheme of the parameters which appear in the LO contribution: 
namely, the top-quark mass. Here, the QCD part of the top-quark mass counterterm is 
defined in the MS scheme. 
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The leading term in the large-m 4 expansion of the two-loop electroweak corrections 
■g(EW) can ^ e f ounc } j n Reference [9], while the hitherto unknown full two-loop result is 

computed in this paper. The sum of and x\ E will only be invariant under an 
electroweak scheme change if it is multiplied by the normalisation factor of the effective 
Hamiltonian, 4aG^/(2\/27rsin 2 9w)- Accordingly, the electroweak renormalisation scheme 
has to be fixed for the parameters in the normalisation factor. 

Since in the electroweak theory not all parameters are independent, we have to specify 
the physical input parameters, and we choose the set 

G F , a, M z , M t , andM H . (2.5) 

Here Gp is the experimental value of the Wilson coefficient relevant for muon decay, a the 
fine structure constant, and Mz the Z-boson pole mass. M t is the top-quark mass, where 
QCD corrections are renormalised in the MS scheme, while the on-shell scheme is used for 
the electroweak corrections. The Higgs mass Mh is essentially a free parameter - its value 
is assumed to be consistent with electroweak precision data. 

For fixed input parameters we can now study the remaining residual higher-order un- 
certainty by using different renormalisation schemes. In the following discussion we will 
make use of three renormalisation schemes: 

• The MS scheme for all parameters, 

• the on-shell scheme for all masses and the MS scheme for all coupling constants, 

• or the on-shell scheme for all masses and the weak mixing angle - the QED coupling 
constant is renormalised in the MS scheme. 

The explicit result for xj: EW ^ is different for each renormalisation scheme. In practise, we 
perform our calculation in the MS scheme and transform our result into the respective 
scheme by a finite renormalisation. 

In all three schemes we renormalise the CKM parameters in the MS scheme and use Gp 
as a normalisation factor for the effective Hamiltonian in Equation (1.1). The parameter Gp 
plays a special role, because it is by itself defined as a Wilson coefficient, of the operator 
Qn = {.^ L lpl^L){(iLl p ^eL) which induces the muon decay in the effective Fermi theory. 
To make this more explicit we introduce the following notation: We denote the Wilson 
coefficient for muon decay by G^ = G^' +G^f W ^ + . . ., where the superscript (0) denotes the 
tree level contribution, (EW) the one-loop electroweak corrections, and the ellipses stand 
for terms beyond second order in the electroweak interactions. By Gp we then denote the 
experimental value of G^ as extracted from muon life-time experiments [11, 12]. If we now 
write the effective Hamiltonian (1.1) in the general form 

a ^ ^ 4Gf a 

' cS= 722^1 

we find 



■Hcff - -7=- . 2 a CvQv - —nTT, • 2 a Xt ^ u ' ( 2-6 ) 

a/2 27rsm 6> w V2 27rsm 6> w 



y<o) _ {EW) _ _ C^Gj EW) (9 7 , 

'-'> <-»> \yfi ) 
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The MS Scheme 



In the MS scheme we use 



9i, 92, v, A, trndyt 



(2.8) 



as fundamental parameters. Here g\ and g 2 are the couplings of the SU(2) and U(l) gauge 
group, respectively, v is the vacuum expectation value of the Higgs field, A the quartic 
Higgs self coupling, and y t the Yukawa coupling of the top quark. All these parameters are 
running parameters, depending on the renormalisation scale //. We fix the initial conditions 
of these parameters by expressing the physical parameter set (2.5) through (2.8) using one- 
loop accuracy 1 and fitting the values of (2.8) to yield the experimental values of (2.5). 

We choose to cancel all tadpole diagrams with a finite counterterm. This results in an 
additional finite renormalisation of all massive quantities - a sample diagram is shown in 
Figure 1. In this way we ensure that intermediate results are gauge parameter independent. 

Masses in the On-Shell Scheme 

As a more well-behaved alternative, we use the on-shell definition of the W^-boson and the 
top-quark mass. Since we performed our calculation in the MS scheme, we have to perform 
a finite mass renormalisation. The necessary renormalisation constants consistent with our 
treatment of tadpole diagrams can be found in [13, 14]. 

In addition, we have to specify the renormalisation scheme for the weak mixing angle. 
We will use the following two schemes: 

• In the on-shell scheme the weak mixing angle is defined by s\y = sin 2 8^~ sheU = 
1 - M^/M§. Here the W-boson mass is calculated including radiative corrections 
from the input parameter set (2.5), which introduces a Higgs-mass dependence. In 
addition, the use of the on-shell value for sin 2 9w implies a finite renormalisation of 
our MS results by including a finite counterterm for sm 2 9\y- It is given in terms of 
the on-shell renormalisation constants for My/ and Mz by 



where the subscript A = implies setting the pole part including the finite subtrac- 
tion, A = 1/e — 7£ + log47r, to zero. The expressions for SM Z and SM^ can again 
be found in [13, 14]. 

• The MS definition of the weak mixing angle, denoted by snd, leads to numerically 
tiny NLO corrections. It is given in terms of s 2 ^ by [15] 




(2.9) 




ee sin 2 6 




1 + 



c 2 w Ana(M z ) 
s 2 s 2 

b W ND 



Ap 



■) 



(2.10) 



For the Higgs boson mass we use the tree-level relation. 
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w<, w s z s w > >w 



t V £ 

Figure 1: Sample penguin, box, and counterterm diagrams. Our tadpole renormalisation results 
in an explicit finite renormalisation of all massive quantities. The right-hand side diagram shows 
a resulting counterterm diagram. 



where a = a MS , and = 1 — s^- The explicit expression for Ap can also be found 
in [15]. 

The numerical discussion of the three different schemes is given in Section 4. 

3 Calculation 

We determine the effective Hamiltonian by computing the relevant Standard Model Green's 
functions in the MS scheme and matching them to the five-flavour effective theory. To this 
end we have to calculate two-loop box and penguin diagrams, samples of which are shown in 
Figure 1. All diagrams reduce to two-loop vacuum diagrams after setting external momenta 
and light masses to zero. The resulting loop integrals are computed using standard methods 
[16, 17]. All this is done in two independent setups: one is using the FeynArts [18] package 
to generate the diagrams and a self written Mathematica program, the other method uses 
a self written Form [19] program. The Feynman gauge £ = 1 is used in both setups. 

The integrals in the effective theory correspond to massless diagrams with vanishing 
external momenta and are exactly zero in dimensional regularisation. The only remaining 
contributions are then products of renormalisation constants and tree-level matrix elements 
of the operators Q u , defined in Equation (2.1), and 



(3.1) 

The evanescent operator E u arises in the context of dimensional regularisation and vanishes 
algebraically in four space-time dimensions. It leads to a non-vanishing finite contribution 
to the Wilson coefficient, proportional to the finite mixing of E u into Q v . The infinite 
operator renormalisation constants are determined from the ultraviolet poles of the matrix 
elements of the operators between external fermion states. They multiply the tree-level 
and one-loop Wilson coefficients of the operators (2.1) and (3.1) and cancel exactly the cor- 
responding spurious infrared divergences of the Standard Model amplitude, thus rendering 
the matching condition finite. 

The use of dimensional regularisation is in general inconsistent with a fully anticom- 
muting 75 matrix in d dimensions, and we use the 't Hooft-Veltman (HV) scheme in our 
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Figure 2: Sample diagrams which imply an off-diagonal field renormalisation. 



calculation. However, problems only arise when computing traces containing at least three 
75 matrices, appearing in the anomalous fermion triangles (see for instance the first dia- 
gram in Figure 1). In all other cases we can use a naive anticommuting 75 (NDR scheme), 
which avoids spurious finite renormalisations required in the HV scheme [20]. 

We have performed our calculation in the MS scheme as described in Section 2. The 
renormalisation of masses and couplings is performed in the usual way. 

In order to ensure the canonical form of the kinetic term for the down-type quarks, 
id^klfidLj, in the effective theory, we perform a finite off-diagonal field renormalisation. 
The exchange of W bosons induces transitions between quarks of different generations (cf. 
Figure 2). We rediagonalise the kinetic term by including a suitable finite part in the 
(matrix-) field renormalisation Z^f-: 

d\7 = Z%d L ^ (3.2) 

where i denotes the generation of the down-type fermion (i = 1, 2, 3). 

The renormalisation leads to a finite result for x[ E . As an additional check we also 
verified that the full result is analytically independent of the renormalisation scale \x. 



4 Numerics 

In this section we present our numerical results and discuss the theoretical uncertainty 
of the branching ratios of the rare Kaon decays. For our numerical analysis we use the 
central values and errors of the input parameters given in Table 1. As discussed in detail in 
Section 2, we use a, Gp, and M z as the basic input parameters for the electroweak theory. 
The mass of the W boson is then not an independent quantity; we calculate its mass using 
the approximate formula given in Reference [21], which includes the state-of-the-art higher 
order corrections. 

Converting the on-shell top-quark mass M t TEV , measured at Tevatron, to the MS scheme 
using three-loop QCD accuracy, we find M t = (mt) = 163.7 GeV. For this con- 

version as well as for the QCD running of M t and a s we use the Mathematica package 
RunDec [22]. 

The electroweak correction term x[ E cancels the scheme and scale dependence of 
the prefactor a/ sin 2 6w up to higher orders in the electroweak interaction. The remaining 
scheme and scale dependence will serve as an estimate of the theoretical uncertainty of our 
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Parameter 


Value 


Ref. 


Parameter 


Value 


Ref. 


My 


91.1876(21) GeV 


[231 


olJMt) 


0.1184(7) 


[231 




155(40) GeV 




&-\M z ) 


127.925(16) 


[23] 


M TEV 


173.3(1.1) GeV 


[24] 


Gp 


1.166 367(5) x 10- 5 GeV -2 


[23] 


m c (m c ) 


1.279(13) GeV 


[25] 


X 


0.2255(7) 


[26] 




0.2315(13) 


[23] 


\v cb \ 


4.06(13) x 10" 2 


[23] 


K+ 


0.5173(25) x 10" 10 


[27] 


p 


n 141+0.029 


[28] 




2.231(13) x lO" 10 


[27] 


fj 


0.343(16) 


[28] 


\ € k\ 


2.228(11) x 10~ 3 


[23] 









Table 1: Input parameters used in our numerical analysis. 



result. To facilitate the discussion, we define the scale and scheme independent quantity 

= aMg) sin 2 e w (n = M z , M H = 155 GeV) 

* a(ji = M z , M H = 155 GeV) sin 2 6w(ji, M H ) ' 

It is formally independent of /j and coincides with X t (fi) at /z = Mz and Mh = 155 GeV. 
We normalise X t to our central value for the Higgs-boson mass, Mh = 155 GeV; as we 
will see below, the dependence on M H is very weak for 115 GeV < M H < 200 GeV. The 
function X t is plotted in Figure 3 for Mh = 155 GeV. Here the dashed line shows the 
LO result. As is clearly visible, the inclusion of the two- loop electroweak corrections (solid 
line) cancels the scale dependence of the electroweak input parameters almost completely, 
up to negligible corrections of 0.02%. 

Next we discuss the dependence of our result on the choice of the renormalisation 
scheme. The difference between the MS and on-shell definition of the parameters sin 2 0\y 
and m 2 , appearing in the LO effective Hamiltonian, amounts to roughly 4% and 7%, 
respectively, leading to a large dependence of the branching ratios on the renormalisation 
scheme, if the two-loop electroweak corrections are not included. In turn, we will see that 
the inclusion of x[ EW ^ cancels this ambiguity almost completely. To get a quantitative 
estimate, we evaluate the function X t numerically in the three renormalisation schemes 
described in Section 2. 

In Figure 4 we show X t in dependence on the Higgs boson mass Mh, where all couplings 
are defined in the MS scheme and all masses in the on-shell scheme. In this scheme the 
NLO electroweak corrections are tiny, of the order of one per mil, even for very large Higgs 
masses. 

In Figure 5 we compare the results in two other schemes. In the left panel we show 
X t , where all parameters are defined in the MS scheme. In the right panel, all parameters 
are defined in the on-shell scheme, apart from a, which is defined in the MS scheme. As 
expected, we observe that for the on-shell definition of sin 2 6w (right panel) the related 
ambiguity is cancelled by a sizeable (~ 4%) two-loop correction, whereas for the full MS 
definition (left panel) the electroweak corrections amount to 1%. 
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1.51 



LO 
NLO 



1.5 



1.49 



1.48 



1.47 



SO 



120 



160 



200 
fi [GeV] 



240 



280 



320 



Figure 3: X t (see text) as a function of //, for Mjj = 155 GeV. The LO result is represented by 
the dashed line, the solid line includes the full two-loop electroweak corrections, which cancel the 
/^-dependence of the LO result almost completely. 

We thus conclude that the on-shell definition of the masses together with the MS 
definition of sin 2 Qw is the best choice of the renormalisation scheme. We can read off 
the maximal difference of the three renormalisation schemes from the two NLO curves in 
Figure 5, right panel - it amounts to 0.27%. For our numerics below, we will take the 
average of the two curves and assign an error of ±0.134% to X t , as resulting from the 
remaining uncertainty of the electroweak correction. In total, using the central values from 
Table 1, we have 



where the first error quantifies the remaining scale uncertainty of the QCD corrections, and 
the second error corresponds to the uncertainty of the electroweak corrections. Here and 
below, we determine the QCD error on X t by varying the scale \i c between 80 GeV and 320 
GeV. Accordingly, our central value of X t is the average of max M X t (/x) and min M X 4 (/i), 
where // G [60 GeV, 320 GeV] . 

Next, let us comment on the validity of the large-m t expansion of the full result, which 
can be gleaned from Figure 5: It is now evident that it is always a bad approximation to 
the full result, as has actually been expected before [9, 29]. 

For convenience we provide an approximate, yet very accurate formula for the NLO 



X t = 1.469 ±0.017 ±0.002, 



(4.2) 
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1.504 



1.496 



1.492 




200 300 

M H [GeV] 



400 



Figure 4: Xt as a function of Mjj. The LO result is represented by the dashed line, the solid line 
shows the result including the full two-loop electroweak corrections. The NLO corrections in the 
limit of large top-quark mass are represented by the dashed-dotted line. 



electroweak correction factor rx = 1 + X\ EW ^ /x[ ^ : 

M t 



r x ^-A + B. C<*»*°»» - D (^) 



where 



A = 1.11508 , 5 = 1.12316, C = 1.15338, D = 0.179454. 



(4.3) 



(4.4) 



It approximates the full result within the limits 160 GeV < M t < 170 GeV to an accuracy 
of better than ±0.05%. 

Finally, we update the theoretical prediction of the branching ratios, including the 
effect of the full two-loop electroweak corrections. After summation over the three neutrino 
flavours the resulting branching ratio for K + — > 7i + uu can be written as 2 [1, 2, 30] 



Br (K + ->• 7T + z/z/(7)) 

= K+Cl + A em) 



(^*)V5£(J>. + «W + 5£* 



. (4.5) 



2 We have omitted a term which arises from the implicit sum over lepton flavours in P c because it 
amounts to only 0.2% of the branching ratio. 
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1.0 
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NLO, s ND ; 
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M H [GeV] 
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1.56 
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1.5 
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LO 
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140 160 
M H [GeV] 



200 



Figure 5: J( as a function of Mh, in two different renormalisation schemes. The dashed lines 
show the LO results, the dashed-dotted lines the LO results including the electroweak corrections 
in the large-mj limit. The full two-loop results are represented by the dotted lines. The left panel 
shows the results where all parameters are defined in the MS scheme. By contrast, in the right 
panel, all parameters apart from a are defined in the on-shell scheme. For comparison, we also 
plot in both panels the NLO result, where all masses are defined on-shell and all couplings in the 
MS scheme. It is represented by the solid lines. 



The parameter 

^>=M^ + H (4e) 

describes the short-distance contribution of the charm quark, where A = \V US \, and has 
been calculated including electroweak corrections, in Reference [6]. The charm quark 
contribution of dimension-eight operators at the charm quark scale \x c [31] combined with 
long distance contributions were calculated in Reference [30] to be 

5P C , U = 0.04 ± 0.02 . (4.7) 

The hadronic matrix element of the low-energy effective Hamiltonian can be extracted 
from the well-measured Ki 3 decays, including isospin breaking and long-distance QED 
radiative corrections [27, 32, 33]. The long-distance contributions are contained in the pa- 
rameters k + , including NLO and partially NNLO corrections in chiral perturbation theory. 
A EM denotes the long distance QED corrections [27]. 

Including the two- loop electroweak corrections to X t , we find for the branching ratio of 
the charged mode 

Br(K + ti + vv) = (7.8lj:^|? ± 0.29) x 10~ n , (4.8) 
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The first error is related to the uncertainties in the input parameters. The main contri- 
butions are (V cb : 56%, p : 21%, m c : 8%, m t : 6%, fj : 4%, a s : 3%, sin 2 W : 1%). The 
second error quantifies the remaining theoretical uncertainty. In detail, the contributions 
are (5P CjU : 46%, X t (QCD) : 24%, P c : 20%, k + : 7%, X t (EW) : 3%), respectively. 

The branching ratio of the CP- violating neutral mode involves the top-quark contri- 
bution only and can be written as 



Br (K L ^«°vv)=K L P^xX . 



(4.9) 



Again, the hadronic matrix element can be extracted from the decays and is now 
parametrised by Kl [27]. There are no more long-distance contributions, which makes this 
decay channel exceptionally clean. 

Whereas the CP-conserving contribution to the branching ratio is completely negligible 
compared to the direct CP- violating contribution within the Standard Model [34], the 
indirect CP-violating contribution is of the order of 1% and should be included at the 
current level of accuracy. This can be achieved by multiplying the branching ratio with 
the factor [35] 

i-^l 1 + p -(W*'>-', (4 . 10) 

V 

where A = V c b/X 2 , and e# describes indirect CP violation in the neutral Kaon system. Tak- 
ing this factor into account, and including again the full two-loop electroweak corrections, 
we find 

Br(K L -> 7T°z/z/) = (2.43+°;^° ± 0.06) x 10~ n . (4.11) 

The first error is again related to the uncertainties in the input parameters. Here main 
contributions are (V c b : 54%, f) : 39%, m t : 6%). The contributions to the second, theoretical 
uncertainty are pf t (QCD) : 73%, k l : 18%, X t (EW) : 8%, 5P C , U : 1%), respectively. All 
errors have been added in quadrature. 



5 Conclusions 

In this paper, we have calculated the complete two-loop electroweak matching corrections 
to X t , the top-quark contribution to the rare decays Kl — > 7T°i/z/, K + — > tt + vv, and B — > 
Xd.sVv. This is in particular important for rare kaon decays: future proposals aim at an 
experimental accuracy of 3% for the branching ratios, while the leading order electroweak 
scheme ambiguity is of similar size. Our calculation reduces the scheme ambiguity in X t 
from ±2% to ±0.134%. The resulting theory uncertainty in the branching ratios is rendered 
from dominant to negligible. 

The absolute corrections are small in a renormalisation scheme where on-shell masses 
and MS coupling constants are used for the electroweak sector. In addition, we analyse 



12 



the convergence in the MS scheme and the on-shell scheme to estimate the remaining 
perturbative uncertainty. 

Our analytic results are summarised by an approximate, but very accurate formula. 
We also give the leading term in a small sin#iy expansion. The full expression can be 
obtained upon request from the authors. 
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A sin 6w Expansion 




-0.02 I ' 

80 120 160 200 240 280 320 

\i [GeV] 

Figure 6: The sinOw expansion of a/Air in the MS scheme as a function of the renor- 

malisation scale fi. The solid line shows the full result, while the dashed line corresponds to the 
leading term of the expansion. 

The explicit expression of the full two-loop electroweak correction X^ EW ^ is too long 
to be given explicitly here. The result significantly simplifies if we expand in the small 
parameter sin 9\y - see Figure 6 for the validity of the expansion. 
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In the MS scheme, after normalising the effective Hamiltonian to Gf, we find 
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X( EW \x t ,a,s NB ,n) = 5>4 + 0(*nd) , (A.l) 



where a = [M H /mf s ) 2 , 

1 1 (1) 

Cl 3a(s t - l)'s t ' C2 (x t - 1)3(^-1)^' 

C3 = 2(rr t -l)3(arr t -l) ¥ ' l( ^' ° 4 = 2{x t - lf{ax t - 1) V?l(i ^ ) ' 

C5 = 2(^-f)3(a^-l) V?l( ^ ) ' C6 = (xt-^axt-l)^ 1 ^ ' 

C7 = 2a^?(x t - 1 l)3(ax t -l) ¥ ' 2( i' ° } ' C§ = a^T l0g2(a) ' 

C9 = 3(s t - l)3(ax t -l) l0g( ^ } ' Cl ° = 2a(x t -l)i t («x t -l) ^ ' 

1 2 f 2 



(x t -l) 2 1Z (a; t -l) 3 



Cl3 = ( Xt -max t -l) l0g(0) ' Cl4 = 2a(x t - lfx t {ax t -l) ^ l0g(a) ' 

C15 = 7 — — r^Li 2 (l - a) , ci 6 = — Li 2 (l-x t ), 

Cl7 = a(x,-l)% t Li2(1 " ^ 



and 



Ai = + (f 6 - 48a)n 2 + (288a - (32 - 88a)7r 2 ):r t + (2003a + 4(4 - 6a - a 2 )7r 2 )x 2 
+ (9a(93 + 28a) - 4a(3 - 2a + 8a 2 )7r 2 )xf 
+ (3a(172 - 49a - 32a 2 ) + 4a(20 - a + 16a 2 )vr 2 )^ 

- (3a(168 + 11a - 24a 2 ) + 4a(45 + 8a 2 )7r 2 )x t 5 
+ 96a7r 2 4 , 

A 2 = - 768x t - (525 - 867a):r 2 + (303 + 318a)x 3 - 195axf , 

A 3 = - 8(95 - 67a + lla 2 )x 2 + 2(662 - 78a - 177a 2 + 40a 3 )x 3 

- (608 + 476a - 595a 2 + 114a 3 )^ + (44 + 188a - 321a 2 + 103a 3 - Sa 4 )x 5 t 

- a(28 - 72a + 33a 2 - 4a 3 )x t 6 , 

A A = + 48 - 10(57 + Aa)x t + 51(29 + 10a)x 2 - (841 + 1265a)x 3 + (308 + M7a)x\ 

- (28 - 40a)4 + Ylax\ , 
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A 5 = + 768 + (816 - 768a)x t + (1240 - 1232a)x 2 - 4(415 + 2a)x 3 t + (311 + 122a)x i t 
+ (145 - 2&Ia)x\ - (36 + 51a)x t 6 + 20ax 7 t , 

A e = + 328x t - (536 + 900a)x 2 t + (208 + 1584a + 670a 2 ):r 3 - a(668 + 1161a + 225a 2 )^ 
+ a 2 (479 + 362a + 28a 2 )x 5 t - a 3 (143 + 42a)x\ + 16a 4 x 7 t , 

A 7 = + 32 - 4(44 - 9a)x t + (384 - 322a - 400a 2 ):r 2 - (400 - 869a - 1126a 2 - 696a 3 )x 3 
+ 2(80 - 488a - 517a 2 - 631a 3 - 264a 4 )x 4 
+ (48 + 394a + 269a 2 + 190a 3 + 882a 4 + 196a 5 ):r 5 

- (64 - 58a - 89a 2 - 95a 3 + 34a 4 + 296a 5 + 32a 6 )^ 

+ (16 - 59a - 79a 2 + 256a 3 - 239a 4 + 57a 5 + 48a 6 )x[ + (1 - a) 3 a 2 (29 + 16a)x 8 t , 
A 8 = + 28a 2 x 2 - 32a 3 x 3 , 

A 9 = - 288 + 36(1 + 8a)x t + 6(647 + 87a)x 2 + 5(55 - 927a - 132a 2 ):r 3 

- (1233 + 98a - 879a 2 - 192a 3 )x 4 + (360 + 1371a - 315a 2 - 264a 3 )x 5 

- 24a(17 - 4a 2 )4 , 

A w = + 32 + 4(-44 + 29a)x t - 12(-32 + 77a + 31a 2 )x 2 

+ 2(-200 + 837a + 767a 2 + 182a 3 ):r 3 - 2(-80 + 625a + 905a 2 + 520a 3 + 82a 4 )x 4 
+ (48 + 1079a + 590a 2 + 1002a 3 + 462a 4 + 32a 5 )x 5 
+ (-64 - 1160a - 501a 2 - 364a 3 - 486a 4 - 72a 5 )^ 
+ (16 + 729a + 1038a 2 + 38a 3 + 238a 4 + 52a 5 )x[ 

- a(192 + 743a + 50a 3 + 12a 4 )^ + 192a 2 x^ , 

A n = + 16x t + 324a; 2 - 36x 4 , 
A 12 = + 216x t - Q72x 2 + 152x 3 , 

A 13 = - 16x t + (16 - 42a)x 2 + (16 + 21a + 60a 2 )x 3 

- (16 - 21a + 45a 2 + 32a 3 )x 4 - a 2 (7 - 24a)x 5 , 

A u = - 32 + (144 - 68a)x t + (-240 + 334a + 332a 2 )a; 2 + (160 - 551a - 660a 2 - 364a 3 ):r 3 
+ a(329 + 451a + 650a 2 + 164a 3 )x 4 + (-48 - a - 59a 2 - 523a 3 - 316a 4 - 32a 5 ):r 5 
+ (16 - 43a - 93a 2 + 255a 3 + 287a 4 + 32a 5 )^ - a 2 (-29 + 42a + 103a 2 + 8a 3 )x 7 t , 

A 15 = - 144(1 - a) 2 x 2 + 144(1 - a) 2 x 3 - 36(1 - a) 2 x\ , 

A 16 = - 32 + 96a + (48 - 32a)x t - 176ax 2 - (16 - 74a)a: 3 + 212ax 4 , 
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A 17 = - 32 + (64 - 100a)x t - 8(4 - 34a - 29a 2 )x 2 - 4a(34 + 170a + 33a 2 )^ 



+ 8a 2 (47 + 51a + 4a 2 K - 16a J (15 + 4a)x t 5 + 32aV 



t ■ 



Here we use 



Li 2 (0 = - 



log(l - &) 



dt. 



and the two-loop functions (pi and y2 2 are given by [17] 



4^^012(2 arcsin( v ^)) , < z < 1 , 
^(21n 2 ^-4Li 2 ^-ln 2 (4^) + §7r 2 ) , z>l, 



and 



(A.2) 



+ 2 In (|(1 + x - y - A)) In (§(1 - x + y - A)) - Inxkiy 
-2U 2 (±(l + x-y-\))-2U 2 (l(l-x + y-\))},\ 2 >0,^+Vy< l i 



^2{x,y) = < 



-A 2 



C\ 2 (2arccos ( ~]^± v ) ) + Cl 2 (2arccos 



l+x-y 
2v^ 



+Cl 2 (2axccos(±=^))|, A 2 <0, v ^+ V ^> 1. 



(A.3) 



Here A 2 = y/l — 1/z and A = a/(1 — x — y) 2 — Axy. The Clausen function is defined by 
Cl 2 (z) = -J°d6\n\2sm{6/2)\. 



B The large-m^ Expansion 

(2) 

The two- loop electroweak corrections to the bbZ vertex, denoted by r b , have been calcu- 
lated in the limit of a large top-quark mass by Barbieri et al. in [36, 37] and were confirmed 
by Fleischer, Tarasov, and Jegerlehner [38] , who found a particularly simple analytic form 
of the results. Buchalla and Buras have extracted from this result the corrections to the 
sdZ vertex, which they used for their analysis of the K — > -nvv decays in [9]. We will now 
take the limit m t — > oo in our complete result for the sdvv transition and compare it with 
the result in [38]. 

Several important points should be mentioned here: As observed in [9], only Z penguin 
diagrams contribute to the sdvv transition in the large-mj limit. The results in [38] have 
been obtained in the so-called "gaugeless limit", where in particular the W boson field 
does not appear. Accordingly, the parameter corresponding to our x t is defined in [38] by 
x t = v / 2G A[ m 2 /(167r 2 ) and will be denoted by x t in our paper. As a consequence, the result 
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for is normalised to G 2 F in [38]. On the other hand, we performed a full Standard 
Model calculation and afterwards took the limit mt — > oo. 

Thus we now take the large-m t expansion of our result, factor out G%, and perform a 
finite renormalisation of the top-quark mass in our LO result, by replacing = M t +SM t , 
in order to transform into the on-shell scheme. Here SM t is given in the large-mt limit by 

„2 



6Mt 



167T 2 S^ 



3 1 1 

- + 1 a 

a 2 



1 
16 



( 4a l/2_ a 3/2^ (a) + A fl 2 lo 



and 
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(B.l) 



9W 



arctanh 



2 - a 



\fifl — A) a 

In this way we reproduce the result in [38] 



+ arctanh 



(B.2) 



_(2),on-shell 



13 , 

9 a- 2a 2 

A 



1 cr 

24 12 



-a tc 



19 



-a H — a In a 



7 o 3 . 

-a a In a 

4 2 



33 

7 — 18a H a 2 



15 



39 



— a H a' 

2 4 



4a 3 Li 2 (l 



a ^i 



It corresponds to the effective Hamiltonian in the limit of large top-quark mass 



H 



AG, 



a 



cff 



-A* 



X t OL 
+ 



■r: 



y/2 27rsin 2 # w '* V 8 ' 4vr32sin 2 ^ w 
Our result in the MS scheme is given by 



(3 + rf))g,. 



T, 



(2),MS 



11 , 

-2 a- 2a' 

A 



1 7 2 

1 a 2 

24 12 



a 
~2 



a + 2a J In a - 

15 39 2 

a H a 

2 4 



-a 2 - ^a 3 ) In 2 a 



4a 3 Li, (1 



7T" 



33 9 

7 — 18a H a 2 

A 



sfag{c 



(B.3) 



(B.4) 



(B.5) 



a (^1 







for /it = M t . It is normalised to GV and thus independent of the tadpole contribution. 
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